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Abstract 

We develop here a simple quantisation formalism that make use 
of Lie algebra properties of the Poisson bracket. When the brackets 
{H, ifi} and {ipi, <fj}, where H is the Hamiltonian and <pi are primary 
and secondary constraints, can be expressed as functions of H and 
cpi themselves, the Poisson bracket defines a Poisson-Lie structure. 
When this algebra has a finite dimension a system of first order partial 
differential equations is established whose solutions are the observables 
of the theory. The method is illustrated with a few examples. 



1. The quantisation of systems with constraints is old as the quantum me- 
chanics itself. The first such problem brilliantly solved was the finding of the 
hydrogen atom spectrum by Pauli in 1926 0. Enforcing the constraints in 
classical mechanics has a satisfactory solution 0, ||, but this is no more true 
in quantum mechanics. The constraints, i.e. a set of functions 

<fi(q,p) = , i = l,2,...,m (1.1) 

restrict the motion of the classical system to a manifold embedded in the 
initial Euclidean phase space and in consequence the canonical quantisation 
rules 

[qi,Pj] = ih6ij 

are no more sufficient for the quantum description of the physical system. 
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In general quantisation is not a well-defined procedure existing today a 
variety of methods which sometimes give different results when applied to 
physical problems, although the starting points are similar from the classi- 
cal point of view. We mention only the people who derive the Schrodinger 
equation by Feynman's path integral method; see for example ||, H ||, who 
find an extra energy term proportional to the Riemann scalar curvature of 
the manifold, even for the simple case of the motion of a particle on the 
n-dimensional sphere. 

The most succesfull method for imposing quantum constraints is that 
found by Dirac M, however nowadays there are some voices who reject it 
claiming that the resulting energy spectrum is incorrect even for simple 
systems || [|. The mechanism found by Dirac was the introduction of a 
new symplectic structure, the Dirac bracket, to handle the second-class con- 
straints and the using of the Legendre multiplicators for finding the true 
Hamiltonian. 

The purpose of this paper is to look at the problem of quantisation with 
constraints from a slightly modified point of view and to show that the new 
proposal leads to correct results. 

When one studies constrained systems one starts with a Hamiltonian, 
H(q,p), and a number of relations of the form (1.1), called primary con- 
straints, which at their turn generate secondary constraints. Let suppose 
that after a finite number of steps the process closes, i.e. no new secondary 
constraints are generated. In the most simple cases one obtains a Poisson 
algebra of the form 

{H(q,p),<p i {q,p)} = Ci<p j (q,p) (1.2) 

{<Pi(q,p),<Pj(q,p)} = Cij(p k (q,p) 

where Cf and are constant structure coefficients. In our opinion this 
Poisson structure is the basic structure for the quantisation procedure. Since 
the Poisson algebra (1.2) transforms by quantisation into a Lie algebra the 
physical observables of the model will be given by the Casimir operators; 
this means that no one of the initial operators transform into a veritable 
observable. 

We applied this idea to the motion of a particle on the n-dimensional 
sphere and we have found that the "Hamiltonian", i.e. the Casimir of the 
corresponding algebra is a quadratic function in the old Hamiltonian and the 
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constraints ||10|| . This quantity is the square of the angular momentum, a 
result which everybody expected to be so. 

We want to extend this method to more general situations than those 
given by Eqs. (1-2) by developing a formalism which makes use of the Lie 
algebra properties of the Poisson bracket. We hope that this formalism will 
solve at least a part of problems encountered in quantisation with constraints. 

More precisely let (ui,u 2 , ■ ■ ■ ,u r ) denote r functions of In independent 
variables (qi, . . . , q n ,Pi, ■ ■ ■ ,Pn) and suppose that all Poisson brackets {ui, Uj} 
can be expressed as functions of (u±, 112, ■ ■ ■ , u r ). In this case these functions 
form a Poisson-Lie structure and any function of (111,112, • • • ,u r ) belongs to 
this algebra. This kind of structure was first introduced by S. Lie who use the 
name of function group [HJ . The full phase space is R 2n with generic point 



(q, p) and the usual Poisson algebra V = (C°°( y R 2r ), {•,•}) is the setting for 
the problem. 

The systems with constraints are good candidates for such structures 
since we start with a Hamiltonian and a number of primary constraints of 
the form (1.1). By taking the Poisson brackets {H,(pi} and {(pi,(pj} they 
generate secondary constraints. Let suppose that this process closes and 
at the end we obtain a finite number of independent dynamical variables 
(ill, • • • , Ur) which describe the dynamics of the constrained system. These 
dynamical variables satisfy a system of equations of the following form 

{Ui,Uj} = fij(ux, . . . ,u r ) (1.3) 

hj = l ; 2,...,r, i<j 

with fij = —fji skew-symmetric functions. If has a power series expansion 
this may have the form 

fij(ui, ...,u r ) = ciij + b k tj u k + cfjUkUi + ... 

In this approach we make no distinction between Hamiltonian, primary 
and secondary constraints, first or second class constraints, all of them are 
simply dynamical variables living in a democratic society, the rules on which 
they obey being the system of equations (1.3). Now, because the Hamiltonian 
is only one of the pairs, we have to solve the problem of integrals of motion 
for a dynamical system governed by Eqs. (1.3). It seems natural to extend 
the classical solution, F is an integral of motion if {F, H} = 0, to the new 
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context by requiring that F is an integral of motion if 



{F, Ul } = 0, z = l,2,...,r 

We remind that the same condition was used by Dirac too [0], but only 
in the new symplectic structure, the Dirac bracket, {-,-}d, an d not in the 
canonical Poisson structure as we do here. Taking into account the Poisson- 
Lie structure defined by Eqs. (1.3) the above equation is equivalent to the 
following system of first order partial differential equations 

i=r dF i=r dF 

1 = 1 1 1=1 1 

j = 1, • • • ,r 

Being a homogeneous system a necessary condition for the existence of a 



non-trivial solution, F ^ ct, is W% 15 



det \{ui,Uj}\ = det \fij(ux, . . .,u r )\ = (1.5) 

The solution(s) of the system (1.4) will depend in general on all dynam- 
ical variables and will play the role of the classical Hamiltonian for non- 
constrained systems, they being the conserved physical quantities of the dy- 
namical system. The classical theory of first order partial differential equa- 
tions tell us that if the rank of the system (1.4) is p then (1.4) may have 
up to n = r — p independent solutions and the easiest way to obtain them 



is by using the characteristic method [[14], [Lq]. The simplest solutions of the 



system (1.4) are called elementary solutions [|14|, |l~5|| , the general solution 
being an arbitrary continuous and derivable function of these elementary so- 
lutions G = G(Fi, F 2 , . . .). By quantisation {•,•} goes into and the 
observables of the theory will be the solutions of the system (1.4). When the 
algebra (1.3) reduces to that of a semi-simple Lie algebra the solutions F^ 
will be the Casimir operators of this algebra and if the respective algebra has 
rank I there will be I Casimir operators by the well-known result of Racah 
|i~6|j . Thus Eqs. (1.3)-(1.4) represent a generalisation of the known powerful 



machinery of representation theory of Lie algebras and give us a method for 
finding the maximal set of commuting observables for a given physical sys- 
tem. Finding the physically relevant operators and their spectra is one of 
the goals of any quantum theory 
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2. In the following we shall illustrate the new method with a few examples 
to show that its content is not void. 

2.1 We consider first the motion of a particle on a n — 1-dimensional 
sphere which is the toy model for testing quantum constrained dynamics 
M, III], O, O . The free Hamiltonian is 



H 



1 V 

2 (p ' p) = y 



where (p,p) denotes the Euclidean scalar product in the n-dimensional space, 
i- e - Si=i Pi • The primary constraint is usually written as 



¥> = {q, q) - R 2 = r 2 



R 2 = 



The Eqs. (1.3) take the form 



{if,H} = 2V {V,H} = 2H 
{<f,V} = 2& + R 2 ) =2r 2 

where V = (q,p) is the secondary constraint. The system of differential 
equations is 

- 2 K§-2(, + tf)f£=0 



OF dF 

2V a^- 2H W = 



(2.1) 



The condition (1.5) is satisfied the dimension of the matrix being odd. By 
applying the characteristic method [14, 15] we have from the last equation 



(p'(t) = <p + R 2 H'(t) 



-H 



The solution is 



p + R 2 = e t 



H = e- 1 

By eliminating t we find that the solution has the form 

F=( V + R 2 )H + g(V) 
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If we use this form in the second equation we get g(V) = —V 2 /2 . Thus 
an elementary solution of the system (2.1) is 

F = UH — V 2 /2 

where U = ip + R 2 = r 2 = (q, q) . 
The " Hamiltonian" will be 

i / n n n \ 

h = uh- v 2 /2 = \ (l> 2 X> 2 - (E^i) 2 ) = 

1 n i n i 

2 Z)(?iPj ~ HPif = 2 E L J = 2 L2 



Thus the quantum observable is the square of the angular momentum [ID 



11, O. Let show that 7i is the good classical Hamiltonian of the problem. 



The Hamilton equations 



have the form 



. _&H . &H 

qj ~Wj Vl ~ ~Wj 

4j =Pj(q,q) -Qj(q,p) 



Pj = -qj {q,v)+Pj (q,p) 

Multiplying the first equation by pj, the second by qj and taking the sum 
we get 

d, x dV n 
Similarly multiplying the first equation by qj we obtain 

Wj = \~^r = _ (q^q)(q^p) = 

which shows that U and V are constant in time and if the constraints are 
fulfilled at the initial time they will be fulfilled at any time. We consider the 
last two relations as a consistency check of the formalism. 

2.2 We consider now a more complicated structure, the functions /y en- 
tering Eqs. (1.4) being quadratic functions. One of the first such a structure 
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is that introduced by Sklyanin in connection with the Yang-Baxter equations 



[17]. The eqs. (1.4) have the following form 

{u 2 , Ui} = hu 3 u 4 {u 3 , Ui} = b 2 u 2 u 4 

{u 4} Mi} = b 3 u 2 u 3 {u 3 , u 2 } = aiUiU 4 

{u 4 , u 2 } = a 2 uiu 3 {w 4 , u 3 } = a 3 u 4 u 2 

a>i and hi , % = 1, 2, 3 being arbitrary complex numbers. The case considered 
by Sklyanin was a\ = —a 2 = a 3 and b\ + b 2 + 63 = 0. The system (1.5) has 
the form 

dF dF OF n 

biu 3 u 4 - h b 2 u 2 u 4 - h a 2 MiM 3 T^— = 

ou 2 ou 3 ou 4 

dF OF OF . , 

-biu 3 u 4 - h a^UA- h b 3 u 2 u 3 —— = (2.2) 

OUi ou 3 ou 4 

dF dF dF 
-b 2 u 2 u 4 - a 3 uiu 4 - h aiuiu 2 —— = 

OUi OU 2 OU4 

dF dF dF 

-b 3 u 2 u 3 - a 2 uiu 3 - aiuiu 2 - — = 

du\ du 2 du 3 

The condition (1.5) is equivalent to 

ai&i - a 2 b 2 + a 3 b 3 = (2.3) 

so in the following we suppose that (2.3) holds. From the first equation we 
have 

u' 2 (t) = b\u 3 u 4 , u' 3 (t) = b 2 u 2 u 4 , u' 4 (t) = b 3 u 2 u 3 

From the first two relations we have that b\u 3 — b 2 u 2 = ct. This suggest to 
look for a solution of the form 

F(ui, u 2 , u 3 , m) = hul - b 2 u\ + g(ui) 

independent of u 4 . The substitution of this F in the second equation (2.2) 
gives g(ui) = a 3 u\ and the first Casimir has the form 

C\ = a 3 u\ — b 2 u\ + b\u\ 

In the same manner one finds the second solution which is 

C 2 = a\u\ - + b 2 u\ 
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2.3 Another quadratic algebra is found in ref. |18| used to describe the 
kinematical symmetry of a spin chain on a one dimensional lattice.lt has the 
form 

{ Ul ,u 2 } = --ul {u 1 ,u 3 } = au 1 (2.4) 

{u 2 ,u 3 } = au 2 {w4,Wi} = z = 1,2,3 

Since w 4 commutes with the other generators a solution of the eqs. (1.4)) is 
of the form /(W4) with / an arbitrary derivable function. The other Casimir 
is 



C = U\U 2 1 — -u\ 



1 

2 

If we perturb the second equation (2.4) to the following form 



{ui,u s } = aui + bu 2 (2.2) 



obtaining a Poisson-Lie structure on the 2-dimensional Galilei algebra [IS 
the Casimir is more complicated and cannot be guessed simply. The charac- 
teristic method gives 



C = auiu 2 1 — b log u 2 — —u 3 



2.4 An other interesting example appears in the construction of Wess- 
Zumino-Witten models on non semi-simple groups |2(| . The algebra has the 
following structure 

{T, J} = {T, Pi} = , i = 1,2 

In general, given a Lie algebra to define a WZW model one needs a bilinear 
form in the generators of the algebra, form which is symmetric, invariant and 
non-degenerate. Usually for semi-simple groups one takes TruiUi, with the 
trace taken in the adjoint representation of the group. For nonsemisimple 
groups this quadratic form is degenerate. By applying our formalism one 
finds easily the two Casimirs 

d = Pi + Pi + 2JT , C 2 = g{T) 
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where g(T) is an arbitrary derivable function of T. Thus the most general 
bilinear form is 

Q = a (P? + P 2 + 2 JT) + bT 2 

where a and 6 are two arbitrary constants, which is the result of Nappi and 
Witten. 

2.5 Now we want to show that finding the spectrum of the hydrogen atom 
is also a problem of quantisation with constraints. The classical Hamiltonian 
is 

2m r 

where m is the reduced mass and k = Ze 2 . H and the angular momentum 
L = r x p are constants of the motion. But these quantities are not enough 
to make the orbit to be closed, and not enough for having a discrete spectrum. 



We quote from Schiff's book |21] 



"The rotational symmetry of H is enough to cause the orbit to lie in some 
plane through O, but is not enough to require the orbit to be closed. A small 
deviation of the potential energy from the Newtonian form V(r) = —(k/t) 
causes the major axis PA of the ellipse to precess slowly, so that the orbit 
is not closed. This suggests that there is some quantity , other than H and 
L, that is a constant of the motion and that can be used to characterise the 
orientation of the major axis in the orbital plane." 

Such a quantity, which we see as a constraint, is the Laplace- Runge-Lenz 
vector. It is proportional with the Di-polar momentum of the orbit and has 
the form 

p X L K 

M = r 

m r 

These constraints generate the first quadratic algebra in quantum physics. 
Indeed we have 

{M u Mj} = --e ijk HL k i,j = 1, 2, 3 
m 

where H and L are the energy and, respectively, the angular momentum. 
The energy commutes with all the other quantities 

{H, Mi} = {H, Li} = 2 = 1,2,3 
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and we have also 



{Mi, Lj} — €ij k M k , {Li, Lj} — €ij k L k , i — 1, 2, 3 

The Eqs. (1.4) have the form 

T dF dF „, OF „, <9F 

-£37; h £27; M 3 - h M 2 - = 

3 dL 2 dL 3 s dM 2 dM 3 

T OF T OF „, OF ^ dF 

L ^r Li ^ 3 +M3 mr Mi dM- 3 -° 

T OF T OF tr dF . r OF n /o rx 

- L ^ + i ^- M ^ + Ml M4=° (2 - 5) 

„ dF ir OF T OF T OF n 

-M 3 - h M 2 - h <xL 3 - aL 2 - = 

dL 2 dL 3 A dM 2 dM 3 

„ dF OF OF 

„ dF tr dF dF T dF 

-M 2 - h Mi- h tt-^iT^ = 

dLi dL 2 2 dM 1 1 dM 2 

where a = 2H / m. Since H commutes with all the other quantities it is in the 
centre of algebra and, such as in the previous example, it will be a Casimir, 
i.e. an observable in the quantum theory. Thus C\ — H — E is a good 
quantum number. 

One can easily see that L 2 and M 2 are separately solutions of the first 
three equations (2.5), but none of them satisfies the last three equations. We 
look for a solution of the form 

F = a L 2 + b M 2 

with a and b some constants. We find from the fourth equation that b/a = 
— 1/a and the second Casimir is 

C 2 = a(L* - «!) 

a 

The third is 

C 3 = L-M 
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If we use the quantum form of M 2 , i.e. 

M 2 = —(L 2 + h 2 ) + K 2 

m 

and take a = 1/4 in the second Casimir, we find the known form of the 
energy levels 

mK 2 

~ ~2h 2 (2c+l) 2 

where c = 0, 1/2, 1, . . . are the eigenvalues of Ci. In conclusion the hydrogen 
atom has a symmetry group, but this is a nonsemisimple one, its Lie alge- 
bra has dimension seven, and, more important, for explaining the discrete 
spectrum is not necessary to invent a dynamical symmetry like 0(4). 

3. The Dirac quantum theory |7j was patterned after the classical theory, 
the " observables" representing constraints must have zero expectation values. 
This requirement is not consistent with the fact that the Poisson brackets 
between Hamiltonian and constraints and between constraints themselves 
may not vanish such as Eqs. (1.3) show. In this paper we have shown that 
this inconsistency disappears if we postulate that the observables are the 
Casimir operators of the algebra (1.3). A consequence of this postulate is 
the following, starting with a Dirac form Hamiltonian 

H D = H + Hiipi 

may be misleading and cause troubles when using it for the description of 
physical systems, the true Hamiltonians being more complicated functions of 
both the old Hamiltonian and the constraints together, as the above examples 
suggest. The lesson to be learnt is that for constrained systems almost no 
one of the initial dynamical variables transforms into an observable. In this 
respect the hydrogen atom is an exception, the reason being that the classical 
Hamiltonian commutes with all the constraints, being in the centre of the 
Poisson-Lie group. 

Acknowledgements. This work was done while the author was a visitor at 
Institute of Theoretical Physics, University of Bern in the frame of the Swiss 
National Science Foundation program " Cooperation in Science and Research 
with Central and Eastern European Countries and New Independent States 
1996-1998. Institutional Partnership". I take this opportunity to thank the 



11 



Swiss National Science Foundation for support. The warm hospitality of 
Professor H. Leutwyler is gratefully acknowledged. 



References 

W. Pauli, Z. Phys. 36 (1926), 336 

E. T. Whittaker, Analytical Dynamics, Cambridge University Press, 
Cambridge, 1937 

J. Marsden and A. Weinstein, Rep.Math.Phys. 5 (1974), 121 
B. S. De Witt, Rev.Mod.Phys 29,(1957),377 
K. S. Cheng, J.Math.Phys. 13 (1972), 1723 
H. Dekker, Physica 103A (1980), 586 

P. A. M. Dirac Lectures on Quantum Mechanics, Yeshiva Univ., N. Y., 
1964 

H. Kleinert and S. V. Shabanov, Phys.Lett. A232 (1997), 327 
J.R.Klauder and S. V. Shabanov, Nucl.Phys. B511 (1998), 713 
P. Dita, Phys.Rev. A56 (1997), 2574 
H. Grundling and C. A. Hurst frep-th/9712052 



E. Abdalla and R. Banerjee |quant-ph /980302T 



S. Lie Math.Ann. 8 (1875), 215. Quoted in ref. [0 at page 322. 

E. Kamke Differentialgleichungen. II, Akademische Verlag, Leipzig 1956 

E. Goursat Legons sur I'integration des equations aux derivees partielles 
du premier ordre, Hermann, Paris 1921 

E.Racah, Rend.Lincei 8 (1950), 108 

E K Sklyanin, Funct.Anal.Appl. 16 (1982), 263 



12 



[18] F. Bonechi, R. Giachetti, E. Sorace and M. Tarlini |q-alg/9610006| 

[19] E. Kowalczyk |q-alg/ 97040091 

[20] C.R.Nappi and E.Witten, Phys.Rev.Lett. 71 (1993), 3751 

[21] L.I.Schiff Quantum Mechanics, McGraw Hill, New York (1955), p. 235 



13 



